function [old_sj_fuel,new_sj_fuel,sj_change_fuel] = change_ms(sj,p,Elast,ind_market,id_fuelcat,group,variation, nmrk)
p_change=(variation.*(group));
%sj_change=(p_change(ind_market==nmrk)'*Marg_p)';
sj_change=sum((p_change(ind_market==nmrk)'.*Elast{:,:,nmrk}),2);
new_sj=sj(ind_market==nmrk).*(1+sj_change);
old_sj_fuel=accumarray(grp2idx(id_fuelcat(ind_market==nmrk)),sj(ind_market==nmrk));
new_sj_fuel=accumarray(grp2idx(id_fuelcat(ind_market==nmrk)),new_sj);

sj_change_fuel=1-(old_sj_fuel./new_sj_fuel);
end